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Abstract 

The accretion-powered pulsar Her X-1 was observed with Suzaku twice in its 
main-on state, on 2005 October 5-6 and 2006 March 29-30, for a net exposure of 
30.5 ks and 34.4 ks, respectively. In the 2005 and 2006 observations, the source was 
detected at an average 10-30 keV intensity of 290 mCrab and 230 mCrab, respectively. 
The intrinsic pulse period was measured on both occasions at 1.23776 s by HXD- 
PIN, after barycentric and binary corrections. The pulse phase-averaged spectra in 
the energy range above 10 keV are well fitted by "Negative and Positive power-law 
times Exponential (NPEX)" model, multiplied by a fundamental cyclotron resonance 
scattering feature at ~36 keV which appears very significantly in the HXD-PIN data. 
The resonance profiles were reproduced successfully by the Lorentzian type scattering 
cross section, rather than by a Gaussian type alternative. The pulse phase-averaged 
HXD-GSO data, covering 50-120 keV, are featureless. However, in a differential 
spectrum between the pulse-decay phase and off-pulse phase, the second harmonic 
cyclotron resonance was detected in the GSO data at ~73 keV, with a depth of 
1.61q;7. This makes Her X-1 a 6th pulsar with established second harmonic resonance. 
Implications of these results are briefly discussed. 
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1. Introduction 

Magnetic field strength is one of tlie important fundamental pliysical parameters of 
neutron stars. Their surface magnetic-field strengths can be most accurately determined by 
measuring quantized electron cyclotron resonances, corresponding to transitions between adja- 
cent Landau levels which are separated by 

E., = n.Q B^2-{l + z)-' keV, (1) 

where B12 is the magnetic field strength in units of 10^^ Gauss, and z is gravitational redshift. 
Since this E^, with B12 ~ 1 falls in the X-ray energy range, accretion-powered X-ray pulsars 
provide an ideal laboratory where we can directly measure Ea-, and hence B12. Indeed, spectral 
absorption features at this resonance, called cyclotron resonance scattering features (CRSFs), 
have so far been detected from more than 15 acceretion-powered X-ray pulsars (e.g., Trumper 
et al. 1978; Wheaton et al. 1978; Clark et al. 1990; Mihara 1995; Makishimal et al. 1990; 
Makishimal et al. 1999; Coburn et al. 2002; Salvo et al. 2004). Using equation (1), the 
surface magnetic field strengths of these pulsars have been found to cluster in a narrow range 
of (1-5) X 10^2 Qauss (Makishimal et al. 1999). 

Some of those pulsars with CRSFs exhibit multiple harmonic absorption features. It 
was reported that 4U 0115+63 has four harmonics (Santangelo et al. 1999; Nakajima et al. 
2006), and X0331+63 has up to the third harmonics (Pottschmidt et al. 2005; Tsygankov et 
al. 2006a; Mowlavial. 2006). In addition, there are objects exhibiting double (fundamental 
and second harmonic) CRSFs, including Vela X-1 (Kreykenbohm et al. 1998; Makishimal et al. 
1999), 4U1907+09 (Cusumano et al. 1998; Makishimal et al. 1999), and A0535+26 (Kendziorra 
et. al. 1994; Grove et al. 1995). 

Since the fundamental and higher harmonic resonances involve somewhat different el- 
ementary processes (e.g., Alexander & Meszaros 1991), measurements of centroid energies, 
depths, and widths of higher harmonics are expected to provide valuable information on the 
physics of electron vs. photon interaction in the accretion column. Nevertheless, we do not have 
sufficient understanding as to, e.g., what controls the relative depths between the fundamental 
and second harmonic CRSFs. Therefore, we still need to enlarge our sample. 

The accretion-powered X-ray pulsar Hercules X-1 (Her X-1) is one of the most studied 
objects of this class, over decades since its discovery in 1972 by Uhuru (Giacconi et al. 1973). 
At ~35 keV, Her X-1 exhibits a fundamental CRSF, which is the first CRSF discovered among 
all pulsars (Trumper et al. 1978). Since then, this ~35 keV CRSF has been studied extensively 
with various X-ray missions, including Ginga (Mihara et al. 1990), BeppoSAX (Dal Fiume et 
al. 1998), RXTE (Gruber et al. 2001, Coburn et al. 2002) and INTEGRAL (Klochkov et al. 
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Table 1. Summary of Suzaku observations of Her X-1. 



Epoch 


Start (UT) 


End(UT) 


Exposure (ks) * 


Position 


1 


2005/10/05 15:12:00 


2005/10/06 10:25:00 


30.5 


XIS nominal 


2 


2006/03/29 18:12:00 


2006/03/30 15:22:00 


34.4 


HXD nominal 



* For the HXD, without dead time correction and ehmination of buffer flush intervals. 



2007). Nevertheless, the presence of the second harmonic CRSF in this object has remained 
controversial. While the observation with BeppoSAX obtained evidence of the second harmonic 
at -Ea = 72 ± 3 keV in the descending edge of the main pulse peak (Salvo et al. 2004), the 
INTEGRAL data did not confirm it in 2005 observations (Klochkov et al. 2007). 

Since the flux of an X-ray pulsar is known to cut off steeply above an energy of ^1.5Ea 
(Makishimal et al. 1999), it is generally not easy to detect the second harmonic CRSF at 
~ 2Ea. The Hard X-ray Detector (HXD; Takahashi et al 2007; Kokubun et al. 2007) onboard 
the Suzaku satelite has realized high sensitivity over a broad energy band, employing Si PIN 
photo-diodes (hereafter HXD-PIN or briefly PIN) and GSO scintillation counts (hereafter HXD- 
GSO or briefly GSO), which cover the 10-70 keV and 50-600 keV energy ranges, respectively. 
When the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) is incorporated, the energy 
range to be covered expands to three orders of magnitude. Thanks to the high sensitivity in this 
broad energy band, Suzaku is expected to settle the issue of the second CRSF in Her X-1. In 
the present paper, we describe detailed pulse phase-averaged and phase-resolved spectroscopy 
of Her X-1 made with the Suzaku HXD, and report on our confirmation of the second harmonic 
CRSF. 

2. Observations 

Her X-1 has three characteristic periods (Nagase 1989); the 1.24 s intrinsic spin period, 
the 1.7 day binary period of the pulsar together with its optical companion HZ Her, and the 
35 day on-off period which is usually attributed to disk precession. In order to observe a high 
flux state of Her X-1, we need to sample so-called main-on phase in the 35 days periodicity, 
and avoid the binary eclipses. 

We observed Her X-1 twice with Suzaku since its launch. The first observation was made 
on 2005 October 5 UT 15:12 through October 6 UT 10:25, and the second on 2006 March 29 
UT 18:12 through March 30 15:22. These dates were both chosen to observe the main-on phase 
and to avoid eclipses, in reference to past observations (Zane et al. 2004; Still et al. 2001). In 
both observations, the HXD was operated in the standard mode, while the XIS employed "1/8 
window" option to improve the time resolution (to 1 s) and to avoid event pile up. 

Figure 1 shows a long-term light curve of Her X-1 obtained by the RXTE ASM. As 
indicated there, the two Suzaku observations both sampled the main-on phase as aimed. 
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Fig. 1. A 2-10 kcV light curve of Her X-1 obtained with the RXTE ASM. The two Suzaku observations 
are indicated with red dashed lines. 



Data reduction 



We analyzed the HXD data prepared via version 1.2 pipeline processing. The data 
screening criteria we employed are as follows: (a) the time after passage through the South 
Atlantic Anomaly should be larger than 500 seconds; (b) the target object should be above the 
earth rim by at least 5°; (c) geomagnetic cutoff rigidity should be greater than 8 GV c~^; and 
(d) the data should be free from "buffer flash" (Kokubun et al. 2007). The screenings yielded 
a net HXD exposure of 30.5 ks and 34.4 ks, in the 2005 and 2006 observations, respectively. 

Although our main objective is to search the HXD data for the second harmonic CRSF, 
we briefly utilize the XIS data as well. Therefore, we retrieved the XIS data of the two obser- 
vations, both processed with version 1.2 pipeline. Figure 2a shows the 0.4-10 keV XIS2 light 
curves. The XIS background, though included in these light curves, is completely negligible 
(~7 X 10~^ c s~^). The XIS data were not acquired in the former half of the 2006 observation, 
due to an operation error. In the 2006 XIS light curve, we find a few occasions of intensity 
decrease. Since the HXD light curves do not show any corresponding feature, they are likely 
to be so-called intensity dips, observed occasionally from Her X-1 (Mihara et al. 1991). 

Panels (b) and (c) of figure 2 are the background-subtracted and deadtime-corrected 
hght curves from HXD-PIN (10-70 keV) and HXD-GSO (50-100 keV), respectively. Thus, both 
observations were free from binary eclipses. The source was detected with an average 10-30 keV 
PIN intensity of 16.9 cnt s~^ and 13.1 cnt s~^ in the first and second observations, respectively. 

We constructed GSO background (bgd_d model) using a method developed by Fukazawa 
et al. (2007) for the 2005 and 2006 observations. We used this GSO background model to drive 
the light curves in figure 2 and to perform standard phase average spectrum analysis in §5.1. 
Although this GSO background model is available only in relatively coarse energy bins, finer 
binnings can be incorporated in the phase-resolved spectroscopy (§5.2), which does not depend 
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Fig. 2. Light curves and period changes of Her X-1 in the 2005 (left) and 2006 (right) observations, (a) 
Background-subtracted 0.4-10 kcV XIS2 hght curve, (b) Background-subtracted and dead-time corrected 
HXD-PIN hght curve in the 10-70 keV energy band, (c) The 50-150 keV HXD-GSO hght curve, obtained 
in the same way as the PIN data, (d) The pulsation period determined by the folding analysis of the PIN 
data at each good time interval, before (black) and after (red) correcting for the binary motion. 

on the GSO background model. To analyze the 2005 data, we utilize PIN background model 
called bgd_a developed by Watanabe et al. (2007), while another model (bgd_d) developed by 
Fukazawa et al. (2007) for the 2006 PIN data. 

4. Timing Analysis 

Since the XIS data have a time resolution of 1 s, we conduct the timing analysis only 
on the HXD data. The arrival time of each HXD event was corrected for the orbital motion 
effect of the Earth around the Sun, and that of the satellite around the Earth, using a Suzaku 
specific tool aebarycen (Terada et al. 2007) and the object coordinates as (a,(5) = (16^57™49?83, 
35°20'32f'6). The bottom panels (black) of figure 2 show the pulse period at each good time 
interval, determined after these corrections by the standard folding analysis of the 10-70 keV 
background-inclusive PIN data. The orbital motion of the pulser is clearly seen. 

As a next step, we corrected the event arrival times for the orbital delay At arising in 
the Her X-1 system, using a formula as 



At 



asm? 



■sm 



27r( 



t 

P orh 



(2) 



Here, a is the semi-mejor axis, i is the inclination, c is the speed of light, t is an event time 
(suzakutime with its origin on 2000 January 1th UT 00:00), Porb is the orbital period, and (po 
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Table 2. Orbital parameters of Her X-1 employed in the present paper (Still et al. 2001). 



Parameter Value 

(a/c) sin i (s) 13.19029 

Porb (days) 1.7001673 

00* 0.064 

" Orbital phase origin, dettned 
using equation (2) at suzaku- 
time 0.0. 

is the phase origin. The values of (a/c) sin i and Porb, as given in table 2, were employed, 
while 00 corresponding to the suzakutime 0.0 was calculated using the phase origin and the 
values of Porb given by Still et al. (2001). When we scanned 0o over a range of 0.0-1.0, just for 
a cross check, periodograms calculated for sufficiently long time intervals (~30 ks) exhibited 
the strongest contrast correctly at 0o ~ 0.064. As presented in figure 2(a) in red, this orbital- 
delay correction has brought all the instantaneous period measurements into a constant value, 
Pq = 1.23776, in both observations, with a typical uncertainty of 1 x 10"^ s. Thus, we quote the 
intrinsic pulse period of Her X-1 as 1.23776 ± 0.00001 s, both on MJD 53646 and MJD 53823. 
Any pulse-period difference between the two epochs is comparable to this error. 

Figure 3 shows energy-sorted and background-inclusive pulse profiles of Her X-1 in the 
two observations. These were obtained by folding the data at the period of Pq. In both 
observations, the pulsation was clearly detected up to 90 keV by the HXD. The measured pulse 
shapes are typical of the main high state (Deeter 1998). Although they are similar between the 
two occasions, a trailing shoulder at phase ~0.2 is more prominent in 2005. 

5. Analysis of the HXD Spectra 

In this section, the HXD spectra obtained in the two observations are analyzed; the XIS 
data analysis will be reported elsewhere. We employ a GSO correction factor, which has been 
introduced by Takahashi et al (2007) to reproduce the Crab spectra by a single power-low model 
in the 70-300 keV energy range. In addition, we introduce 1% systematic errors in all spectral 
fitting analyses, to refiect typical uncertainties in the current instrumental calibration. The 
model normalization is constrained to be the same between HXD-PIN and HXD-GSO, while 
allowed to take different values between the HXS and HXD. This is because the 1/8 window 
option of the XIS introduces some uncertainties in the absolute source flux. 

5.1. Phase-averaged spectra 

Figure 4 shows the background-subtracted 0.1-100 keV spectra of Her X-1, obtained by 
the XIS, HXD-PIN and HXD-GSO. Since the 2005 and 2006 data give very similar spectra, 
we co-added them together in flgure 4. The background models (described in §3) used for 
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Fig. 3. Energy-sorted and background-inclusive pulse shapes of Her X-1 obtained in the 2005 (left) and 
2006 (right) observations. The data were folded by the barycenter- and binary-corrected pulse period, 
Po =1.23776 s. The phase zero (0 = 0) is adjusted to the main central peak in the PIN energy band. 

the PIN and GSO data are also shown. From the XIS data, we subtracted the night earth 
backgrounds. After subtracting the background, the summed data yield an average 1-50 keV 
flux of 6.0 X 10~^ erg cm~^ s"^. This value was derived by fitting the XIS, HXD-PIN, and HXD- 
GSO spectra simultaneously by an empirical model consisting of an NPEX continuum (to be 
described later), and a low-energy blackbody to reproduce the soft X-ray excess. Assuming 
an isotropic emission, this implies a 1-50 keV luminosity of 3.1 x 10'^^ erg s~^ at a distance of 
D = 6.6 ±0.4 kpc (Reynolds et al. 1997). This luminosity is typical of Her X-1 in the main 
on state. For reference, the first and second observations yielded the 1-50 keV luminosity of 
3.6 X 10^^ erg s~^ and 2.5 x 10^"^ erg s~^, respectively. 

The HXD spectra in figure 4 are expanded in figure 5a, and are normalized in figure 
5b to the Crab spectra acquired with the HXD on 2005 September 15 at the HXD nominal 
position. Thus, the source intensity averaged over the two observations was ~250 mCrab in 
the 10-30 keV band. Independent analyses of the 2005 and 2006 data gave the same quantity 
as ~290 mCrab and ~230 mCrab, respectively. In the Crab ratio, we observe a broad dip at 
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Fig. 4. Background-subtracted and pulse phase-averaged spectra of Her X-1, obtained by the XIS, 
HXD-PIN, and HXD-GSO. Spectra from the 2005 and 2006 observations are summed up. The non-Xray 
background models of PIN and GSO are also shown. 

~36 keV, to be identified later with the fundamental CRSF. 

To represent the spectral continuum, we employed so-called Negative and Positive power- 
law times Exponential (NPEX) model, expressed as (Mihara 1995; Makishimal et al. 1999) 



Here, E is the energy, f{E) is the photon number spectrum, Ecut is a cutoff energy, A and B 
are normalization factors, while cti > and a2 > are two photon indices. We adopt 02 = 2 
for the positive power-law component, so that the second term represents a Wien hump in a 
saturated inverse Compton spectrum. First, we fitted the data by this NPEX continuum, but 
the fit was not acceptable with a reduced chi-square of ~11. The residuals, shown in figure 5c, 
exhibits a strong dip around 36 keV. This feature, already noted in figure 5b, is interpreted as 
the CRSF established through the past studies (§1). 

As a next step, we multiplied the NPEX continuum by a factor e~'^, where S represents 
the cyclotron scattering cross section given as 



with Ea., D, and W being the energy, depth, and width of the resonance (e.g., Clark et al. 
1990, Makishimal et al. 1999). This model has given a fully acceptable fit with a reduced 
chi-square of 0.51. This value is apparently too small, suggesting a slight over-estimation of the 
systematic error, but this simply makes our subsequent analysis more conservative. Histograms 
in figure 5a display this best fit model, and figure 5d shows residuals between the data and the 
model. The obtained resonance energy = 35.9l[JJ keV, with its depth D = 1.2l[J;} and width 
W = 12.2~^\'l, is generally consistent with the past measurements. The best-fit parameters are 
summarized in table 3. 



f{E) = {AE-'^^+BE^'^') X exp(-E/E,,0- 



(3) 




(4) 
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Fig. 5. (a) An expanded view of the HXD spectra presented in figure 4. Histograms in the figure represent 
the best fit model consisting of an NPEX continuum and a fundamental cyclotron resonance, (b) The same 
spectra, divided by those of the Crab Nebula, (c) Residuals from an NPEX fit. (d) Residuals between the 
data and the best fit model in panel (a), (e) Those when a second harmonic CRSF is included. 

In the same way, we analyzed the 2005 and 2006 spectra separately, and obtained the 
results as presented in table 3. Thus, the NPEX times single CRSF model is successful on the 
individual 2005 and 2006 data as well. The value of Eg^ is consistent, within errors, between 
the two observations, in which the 10-30 keV counts was different by ~30±10%. 

In figure 5b, the Crab ratio suggests a shallow structure at ~70 keV suggestive of the 
second harmonic CRSF (§1). In order to quantitatively examine this possibility, we introduced 
a model consisting of an NPEX continuum, multiplied by two CRSF factors both of the form 
of equation (4). The energy i?a2 and width W2 of the second CRSF factor were fixed at 2£'ai 
and 2Wai, respectively (with the suffix 1 specifying the parameters of the first CRSF), while 
its depth D2 was left free to vary. However, as presented in figure 5e, this new model did 
not significantly improve the fit to the 2005+2006 spectrum: the fit chi-square decreased only 
by 1.2, while the degree of freedom changed from 69 to 68. Therefore, the second harmonic 
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resonance is insignificant, with the 90% upper hmit being D2 < 1.5. This is not surprising, 
since figure 5d reveal httle evidence for a negative feature at ~ 2E.^ = 72 keV. 

We repeated the same analysis, namely the NPEX times double CRSF fit, on the 2005 
and 2006 data separately. The results were essentially the same as before; the chi-square 
changed from 45.4 (z/ = 69) to 44.3 (z/ = 68) for the 2005 data, and from 58.4 {u = 69) to 58.1 
(z/ = 68) for the 2006 data, implying an insignificant improvement in both cases. 

5.2. Phase-resolved spectra 

In many pulsars, the continuum spectra, as well as the cyclotron resonance energies 
and depths, are known to depend significantly on the pulse phase (e.g., Klochkov et al. 2007). 
Especially in the phase resolved spectra of Her X-1, Salvo et al. (2004) reported a second 
CRSF at the descending edge of the main pulse peak (§1). We therefore proceed to pulse 
phase-resolved spectroscopy, using the 2005 and 2006 data summed up to increase statistics. 
In the present case, the phase-resolved analysis has two additional merits: it allows us to avoid 
uncertainties in the HXD background models because we can take direct spectral differences 
between different pulse phases, and to use finer bindings of the GSO data than is specified by 
the current GSO background model (§2). 

Referring to figure 3, we adopted five on-pulse phases as 0.8-0.9, 0.9-1.0, 1.0-1.1, and 
1.1-1.2, while the off-pulse phase as 0.2-0.8. The phases 0.8-0.9 and 1.1-1.2 correspond to 
the soft leading and descending shoulders of the main peak, while those of 0.9-1.0 and 1.0-1.1 
correspond to the leading and descending halves of the main hard peak, respectively. The phase 
resolved spectra, obtained by subtracting that of the off-pulse phase, are shown in figure 6. In 
these spectra, the PIN and GSO backgrounds are considered to cancel out with a considerably 
higher accuracy than is available via the model background subtraction (~5% for PIN and ~2% 
for GSO), because the pulse period of 1.24 s is sufficiently shorter than the time scales of typical 
background variations (minutes to hours), and because Her X-1 is not bright to cause significant 
(> 1%) increases in the HXD dead times. The bottom half of each panel in figure 6 presents 
the data, divided by the best-fit NPEX times CRSF model determined by the phase-averaged 
spectra; the model is meant to provide a rough standard for the phase-resolved spectroscopy. 
Thus, some phase-resolved specra are approximately represented by the phase-averaged NPEX 
times CRSF model, but in general the spectral shapes depends significantly on the phase. 

We fitted the four pulsed-component spectra by an NPEX continuum and a fundamental 
CRSF model, which has given a satisfactory fit to the phase-averaged data. Table 4 summarizes 
the best fit parameters for individual on-pulse phases, and figure 7 shows the fit results for the 
particular phase = 1.0 — 1.1. In the case of 0=0.8-0.9, 0=0.9-1.0, and 0=1.1-1.2, we obtained 
acceptable fits by this model, with xl = 1.05, 0.51 and 0.62, respectively. The fundamental 
CRSF is highly significant in all these phases, and the derived resonance energy, Eg^ = 40.6 keV 
(0 = 0.8 - 0.9), = 39.0 keV (0 = 0.9 - 1.0), E^ = 36.5 keV (0 = 1.0 - 1.1; see figure 7c), 
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Fig. 6. Phase-resolved spectra obtained by subtracting the ofF-pulse (phase 0.2-0.8) spectrum from 
on-pulse ones; (a) =0.8-0.9; (b) 0=0.9-1.0; (c) </> =1.0-1.1; (d) </> =1.1-1.2. The 2005 and 2006 ob- 
servations are summed up. The bottom half of each panel shows ratios of the data, to the best-fit NPEX 
times single CRSF model obtained by the phase-averaged spectra. 

and 34.8 keV (0 = 1.1 — 1.2), stays within ~ ±10% of that determined with the phase-averaged 
2005-1-2006 spectra. However, as presented in figure 7c, the model failed to reproduce the 
= 1.0 — 1.1 spectra {xl = 1-67), due to significant negative residuals at ~70 keV in the GSO 
data which is visible even in the ratio between the data and the phase-averaged NPEX times 
single CRSF model (figure 6c). This can be interpreted as the second harmonic resonance. 

We introduced a second CRSF factor, as already performed in §5.1, again with the 
second resonance energy £'^2 fixed at 2£'a. As shown in table 4, the second CRSF factor hardly 
improved the fits to the (p = 0.8 — 0.9, (p = 0.9 — 1.0, and </> = 1.1 — 1.2 spectra. However, the 
fit to the (p = 1.0 — 1.1 spectra has been drastically improved from xt = 1-67 to xt = 1-07, and 
has become acceptable. An F-test indicates that this improvement being cause by chance is 
~ 1.9 X 10^^. The second harmonic depth is obtainded as D2 = 2.4tij (figure 7d). We therefore 
conclude that the second harmonic CRSF is significantly present in the 2005-1-2006 spectrum 
of the pulse phase = 1.0 — 1.1. 

In order to examine how the present data can constrain the second resonance energy 
£'a2; we repeated the fitting to the = 1.0 — 1.1 spectra incorporating the two CRSF factors, 
but scanning i?a2 over 55-95 keV. We fixed W2 at 2Wi. The behavior of the minimum x^, 
achieved at a given value of -E'a2, is presented in figure 8. Indeed, the chi-square has reached 
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Fig. 7. (a) Pulsed-component spectra of Her X-1, obtained by subtracting the data of the off-pulse phase 
((/) = 0.2 — 0.8) from those of the declining phase {4> = 0.0 — 0.1) of the main pulse peak. Data from the 2005 
and 2006 observations are summed up. The histograms in the figure are the best fit model consisting of 
an NPEX continuum and double cyclotron resonances; corresponding residuals arc displayed in panel (c). 
(b) Residuals when only an NPEX continuum is fitted, (c) Fit residuals from the NPEXxCRSF model, 
with Ea = 36.5 keV. (d) Those with an NPEX and two CRSFs, with Ea2 fixed at 2Ea and W2 fixed at 
2Wi. (e) The same as panel (d), but Ea2 is set free. 

minimum at ~70 keV, where the second CRSF is expected to appear. This ensures that the 
GSO spectrum at this pulse phase has a significant negative feature at an energy that is close 
to 2i?a = 73.0 ±1.0 keV at the same pulse phase. This NPEXx2CRSFs model yields the second 
harmonic resonance energy as E^2 = 70.2tli, and its depth D2 = 1.6^1]:? (figure 7e). 

So far, the second resonance width W2 was fixed to twice that of the fundamental. When 
W2 is made free to vary, we obtained W2 = 29.3l^ 9 by constraining Ea.2 = 2£'a, or W2 = 21.9l^ g 
by allowing E^2 to vary freely. Although the upper limit on W2 thus becomes unbound, the 
90% error range of W2 still includes 2Wi, confirming the consistency of our assumption. 
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Fig. 8. Decrement in the fit chi-square compared to the NPEXxCRSF fit. achieved by adding a second 
CRSF factor at various energies. The same spectra as presented in figure 7 are used. The width W2 is 
fixed at 2Wi. The chi-square of 83.5 obtained by an NPEX and a fundamental CRSF model is also shown. 

6. Discussion 

6.1. Pulse periods 

We observed Her X-1 twice with Suzaku, in 2005 October and 2006 March, and obtained 
high-quality data with the HXD and XIS over an extremely broad energy band (0.1-100 keV). 
On both occasions, Her X-1 showed X-ray intensity and spectra typical in the main-on state of 
its 35-d cycle. 

The intrinsic pulse period was determined with HXD-PIN at 1.23776 ±0.00001 s on both 
occasions. This value agrees with the pulse period history of Her X-1 reported by Staubert et 
al. (2006). 

6.2. The fundamental cyclotron resonance 

In the pulse phase-averaged HXD-PIN spectra obtained on both occasions, the funda- 
mental CRSF has been detected clearly at ~ 36 keV. Using the NPEXxCRSF spectral model, 
we have quantified the resonance parameters (table 3, table 4). 

Based on multiple observations of Her X-1 with RXTE and other missions, Gruber et 
al. (2001) and Staubert et al. (2007) argued that its fundamental CRSF energy exhibited long- 
term variations, possibly related to luminosity changes: after having stayed at ~ 35 keV until 
1991, it increased after 1991 to nearly 40 keV, and then it appears to be gradually returning 
to the previous level. The value of E^, = 35.9 ± 0.3 keV, which we obtained using the pulse 
phase-averaged 2005+2000 data, is consistent with the pre-1991 value compiled by Staubert et 
al. (2007), suggesting that the CRSF energy-change episode, whatever the cause be, has come 
to an end. 

Because HXD-PIN has a higher sensitivity at ~36 keV than proportional counters, 
and better energy resolution (~ 3 keV in FWHM) than inorganic scintillator instruments, the 
present results provide one of the best-quality data on the fundamental CRSF of Her X-1. 
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Importantly, the 10-50 keV PIN spectra bearing the CRSF, including both the phase-averaged 
and phase- resolved (0 = 1.0 — 1.1) ones, can be reproduced successfully by combining the NPEX 
continuum of equation (3) with the scattering cross section of equation (4); the CRSF profile 
has been resolved with a reasonable accuracy. We therefore reconfirm the appropriateness of 
this modeling, which has been used in many of the past studies of CRSFs including those with 
Ginga (Clark et al. 1990; Mihara et al. 1990; Mihara 1995; Makishimal et al. 1999), RXTE 
(Nakajima et al. 2006), INTEGRAL (Tsygankov et al. 2006a), and Suzaku (Terada et al. 2006). 

As an exercise, we replaced the cyclotron absorption factor exp(— 5) {S referring to 
equation 4) with a Gaussian absorption factor used in some other studies (e.g, Staubert et 
al. 2007; Klochkov et al. 2007); exp[—aexTp{{E — EcY /2a^}], where a, Ec, and a are free 
parameters. When this factor is applied to the NPEX continuum, the fit to the phase-averaged 
HXD (2005+2006) spectra worsened to x^ = 0.70 (i/ = 69), from that obtained using equation (4) 
ixl = 0.52 for u = 69; table 3). Although this fit is still acceptable (due possibly to an over- 
estimated systematic error), the Gaussian absorption has failed to reproduce the </> = 1.0 — 1.1 
spectra with xt ~ 2-72 (z/ = 48), even when two of them (for the two harmonic CRSFs) were 
incorporated. Therefore, at least for Her X-1, the Lorentzian-like form of equation (4), which 
is based on the classical cross section of cyclotron resonance (Clark et al. 1990), is considered 
more appropriate than the alternative Gaussian cross section. 

The above argument may need an important remark: the very deep fundamental CRSF 
of the transient pulsar X0331+53 (V0332-I-53), at about 28 keV, cannot be described adequately 
if using a single form of equation (4) (Makishimal et al. 1990; Nakajima et al. 2006; Tsygankov 
et al. 2006a). Instead, a nested pair of such absorption factors with different widths (e.g., in 
terms of Gaussians; Pottschmidt et al. 2005) may be needed. 

6. 3. Possible origins of the resonance width 

What is the origin of the relatively large width of CRSFs, which is generally expressed 

as 

= (0.2-0.5)^a (5) 

(Makishimal et al. 1999; Nakajima et al. 2006)? It has often been argued that W is determined 
mainly by thermal Doppler effects in the accretion column, including those associated with 
electron motion along the magnetic field lines (e.g., Dal Fiume et al. 1998; Cusumano et al. 
1998; Coburn et al. 2002). However, the failure of the Gaussian optical depth to reproduce 
the present Suzaku data casts doubt on this interpretation. Furthermore, as pointed out by 
Makishimal et al. (1999), the values of W of various X-ray pulsars often exceed the expected 
thermal broadening. According to Nakajima et al. (2006), Wi of the transient pulsar XOl 15+63 
increased by a factor of ~ 5 as its luminosity changed from 2 x 10'^^ erg s~^ to 5 x 10^"^ erg s~^, 
but the temperature of the emission region, oc Ecut of equation (3), increased by only ~ 20%. A 
similar luminosity dependence of W may be visible even between the present two observations. 
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These observed changes in W also argue against the thermal broadening scenario. 

In explaining W, an obvious alternative is pulse phase- dependent changes in E^,. 
However, in our table 3 and table 4, Wi is not necessarily smaller in phase-resolved spectra 
(except in = 1.1 — 1.2) than in the phase-averaged ones. Furthermore, the superposition of 
narrow features with different over the pulse phase would not produce wide Lorentzian-like 
wings. Therefore, this interpretation is not likely, either. 

In classical electrodynamics, the width W in equation (4) represents damping effects on 
the gyrating electrons. In quantum mechanics, the resonance width in this formula generally 
reflects the finiteness of the life time of excited states, as TV ~ hA, through uncertainty principle: 
here, A is transition rate of electrons out of the relevant excited state, and h is the Planck 
constant. This A, in turn, is determined by three competing de-excitation processes; namely, 
coUisional, spontaneous radiative, and stimulated (induced). Among them, the spontaneous 
transition rate Arad (so-called Einstein's A-coefficient) is given as 

A,ad ~ 3.8 X lO^'Bu 

(Meszaros 1992), so that the resonance width due to spontaneous emission, namely natural 
width, becomes 

W^nat ~ hK^d = l5Bu eV = 1.3 X 10^=^^^ (6) 

where the final form employs equation (1) neglecting the gravitational redshift z. This is still 
inadequate to explain equation (5). A brief calculation shows that the coUisional de-excitation 
is even less effective. 

Let us finally consider the stimulated emission, of which the rate is given as Agt = KtJa, 
where Ast is so-called Einstein's B coefficient, and Ja is the radiation energy flux per unit photon 
frequency at the resonance energy E^^; this effect was considered by Alexander & Meszaros 
(1991). As estimated very crudely in Appendix, Agt in Her X-1 could be some two orders of 
magnitude larger than A^ad- Then, from equation (6), the resonance width due to stimulated 
emission, Wst = VFnat(Ast/Arad)5 could be a considerable fraction of E^, and might potentially 
provide an explanation to the observed width. An advantage of this picture is that we can 
naturally explain the observed positive dependence of W on the X-ray luminosity (Nakajima 
et al. 2006), because we expect W^t oc Agt oc Ja- On the other hand, a caveats is that equation 
(7) in Appendix predicts the W/Eg, ratio to decrease steeply toward the higher-field objects, 
although such a scaling is not necessarily supported by observation. Furthermore, it is not 
obvious if the absorption feature can be really produced in such a condition as Ast/Arad ^ 1, 
which would usually emhance emission. We leave this intriguing issue to further studies. 

6.4- The second harmonic resonance 

The pulse phase-averaged HXD spectra did not require the second harmonic resonance, 
beyond a 90% upper limit of D2 < 1.5. However, in the pulse phase-resolved spectra which 
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Fig. 9. The best-fit vfv model spectra. Tlie black line shows the best-fit NPEX continuum multiplied by 
a fundamental CRSF determined using the 2005+2006 pulse phase-averaged spectrum. The red line shows 
the best-fit NPEX continuum multiplied by double CRSFs specified by a 2005-1-2006 pulse phase-resolved 
spectrum (0 = 1.0 — 1.1). 

cover the descending edge of the main pulse (0 = 1.0 — 1.1), we have successfully detected the 
second CRSF with D2 = ^Ati'l (under the constraints of Ea.2 = 2-E'a and W2 = 21^2) • Its energy, 
when allowed to vary, becomes i?a2 = yielding Ei^/Ea, = 1.9 ±0.2 which satisfies the 

1:2 harmonic ratio. Here, the value of E^^ can be taken either as the phase-averaged value, or 
that at this particular pulse phase, because they agree within 2%. 

As evidenced by the large chi-square decrements (figure 8), the present detection of the 
second harmonic feature is statistically highly significant. The feature cannot be artifacts due 
to systematic errors in the background subtraction (§5.2), since background systematics can be 
neglected in these "on-pulse minus off-pulse" spectra. Another concern is the uncertainty in 
the GSO response, toward lower energies (Kokubun et al. 2007). We confirmed that the second 
harmonic resonance remains significant even when we discard the GSO energy ranges below 60 
keV, although errors on the second resonance parameters increase. We therefore conclude that 
the detected feature is real. 

Evidently, HXD-GSO played a major roll in detecting the second CRSF, which lies out- 
side the HXD-PIN energy range. Nevertheless, the PIN data also contributed to its detection, 
because the inclusion of the second resonance factor has considerably decreased fit residuals in 
the PIN energy range as well (figure 7c, d, e). In figure 9, we show the inferred best-fit model 
spectrum for phase-averaged with a fundamental CRSF and for phase-resolved with double 
CRSFs. 

The second CRSF of Her X-1 was first reported by BeppoSAX, from a main-on state 
observation in October 2000 (Salvo et al. 2004). It appeared in the descending edge of the 
main pulse peak, with the resonance energy at £'a2 ~ 72 ± 3 keV and its width W2 ~ 11 ± 1 keV. 
Our results agree with those from BeppoSAX, with respect to the resonance energy, as well as 
the particular pulse phase where it becomes significant. These results hence make Her X-1 a 
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Fig. 10. (top) Fundamental CRSF energies presented as a function of the pulse phase, obtained by 
subtracting the interpulse spectrum from on-pulse ones using the 2005+2006 data. The 10-20 keV pulse 
profile (also summed over 2005 and 2006) is superposed, (bottom) Fundamental CRSF depths presented 
in the same manner. 

6th pulsar with the second harmonic resonance. The negative detection of this feature with 
INTEGRAL, in another main-on state observation conducted in July-August 2005 (Klochkov 
et al. 2007), is probably due to insufficient statistics. 

The present results suggest that the second harmonic resonance is a rather common 
feature of accreting X-ray pulsars, particularly when they are luminous (Nakajima et al. 2006). 
Even though the cross section of photon vs. electron interaction in a strong magnetic field 
is much larger at E^, than at 2£'a, the two-photon effect by Alexander & Meszaros (1991) 
ensures that the second harmonic feature can appear as strong as the fundamental. That is, 
the fundamental resonance acts essentially as a photon scattering process, due to very short 
lifetime of the excited state. In contrast, the second resonance is expected to act as a pure 
absorption, because an electron excited by two Landau levels will return to the ground state by 
emitting two photons with energy ~ Eg,; these photons fill the fundamental feature, and make 
it shallower. 

6. 5. Pulse-phase dependence of the resonances 

Figure 10 shows the measured E^i and Di as a function of the pulse phase; numerical 
values are given in table 4. The results generally agree with the INTEGRAL measurements 
(Klochkov et al. 2007). Thus, in the particular pulse phase (0 = 1.0 — 1.1) where the second 
CRSF was detected, the fundamental resonance is also very deep showing a large value of Di, 
together with a relatively low E'ai (though not the lowest). Then, how this particular pulse 
phase relate to the rotational phase of the pulsar? 

Although the rotational modulation of E'ai does not allow straightforward explanation, 
we can think of a simple physical effect related to it; the luminosity-dependent change in E^i. 
In the transient pulsars X0115+63 (Mihara, Makishima & Nagase 2004; Nakajima et al. 2006) 
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and X0331+53 (Tsygankov et al. 2006b; Nakajima 2006), -Eai "was found to decrease as the 
luminosity becomes high enough (several times 10'^^ erg s~^ or higher). This is presumably 
because the accretion column then gets taller, so the resonance photosphere, as seen from end- 
on directions, appears at a higher altitude in the column where the dipole field intensity is 
lower (Mihara, Makishima & Nagase 2004). If this effect is also taking place in Her X-1, we 
expect to observe a relatively low value of E^, when looking onto one of the accretion poles. In 
contrast, when viewing an accretion column relatively side on, we will sample various heights 
along it, and to measure a higher value of E^. 

Another candidate mechanism to produce the phase-dependent change in Ea is the angle- 
dependent relativistic effect pointed out by Harding and Daugherty (1991). For Ea ~36 keV, 
this effect makes the exact value of Ea to decrease by ~3% at the pole-on phase than at the 
side-on phase. Although this shift is small, it works in the same sense as the previous one. 

In addition to the above two mechanism, yet another effect may also produce the phase 
dependence of E^^; namely, bulk Doppler effect in the accretion column. Since the accretion 
flow onto a pulsar has a typical velocity of Vfi ~ c/3, we expect the post-shock plasma in the 
radiating accretion column to have a bulk flow of ~ ^^£1/4 ~ 0.1c, at least just beneath the 
standing shock surface. Then, the longitudinal Doppler shift associated with this bulk flow will 
make i?a lower by ~ 10% when we look onto the column, while the effect will vanish at side-on 
phases. Actually, Terada, Ishida & Makishima (2004) invoked such Doppler shifts in accretion 
columns of magnetized white dwarfs, and successfully explained the unusually strong atomic 
emission lines observed from several objects of that kind. In X-ray pulsars, this mechanism is 
expected to work in the same sense as the above two effects, thus enhancing one another. 

From these considerations, we tentatively conclude that the descending half of the main 
pulse, where Eg, is relatively low, corresponds to the phase where we are looking end-on into one 
of the two accretion poles. This agree with theoretical pulse-profile decomposition by Leahy 
(2004). According to this work, the descending half corresponds to the pole-on phase of one pole; 
the pulse profile becomes asymmetric with respect to this phase, due to an addition of fan-beam 
emission from the other pole which reaches us after affected by gravitational light bending. If 
this phase assignment is correct, we are to conclude that the second harmonic feature becomes 
most prominent when we are observing one pole from and end-on aspect. Further comparison 
with other pulsars, including in particular X0115+63 and X0331+53, would be of great value, 
though beyond the scope of the present paper. 

7. Conclusion 

In this paper, we analyzed the HXD data from two Suzaku observations of Her X-1. The 
fundamental resonance was clearly detected at ~36 keV in the pulse phase-averaged and phase- 
resolved spectra. The Lorentzian-like form of the resonance was found to be appropriate than 
the alternative Gaussian cross section. The second resonance feature, though absent {D2 < 1.5) 
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in the phase-averaged spectra, was detected at ~73 keV in the pulse phase-resolved spectra at 
the descending edge of the main peak. 

The authors are grateful to all the members of the Suzaku Science Working group, for 
their help in spacecraft operation, instrumental calibration, and data processing. 

Appendix: Estimation of induced emission 

Assuming a local thermal equilibrium, detailed balance relates the Einstein's A and B 
coefficients as Ast = Ai.ad(^c)^/(2i?f), where the quantities are defined in subsection 6.3. Then, 
the stimulated transition rate, relative to the spontaneous one, can be written as 

Ast/Arad=(M'^a/2^f . (7) 

Let us next estimate the spectral fiux density Ja at 36 keV. The observed fiux density 
per unit energy at E^, before the CRSF factor is applied, is derived from our best fit model to 
the phase-resolved spectra as 3.6 x 10~^° erg s~^ cm~^ keV~^. Multiplying this by AnD^, where 
D = 6.6 kpc (subsection 5.1) is the distance, we obtain a monochromatic luminosity aX E = E^^ 
as / = 1.8 X 10'^^ erg^^ s"^ keV^^. Further assuming that the two polar emission regions of Her 
X-1 can be approximated each by a hemisphere of radius r ~ 100 m, we obtain 

47rr^ 

= 3.9 X 10^0 ( ( -^^) \eV s-i cm-^ Rz"' 
VlOOm/ V36keV/ 

after transforming units from erg keV~^ to keV Hz~^. This finally yields 
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Table 3. Best-fit parameters of the phase averaged spectra of Her X-1. 
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JNP represents JNFEX continum model, (Jl represents a toundamental CKbit'. (J12 represents ioundamental and 
second CRSFs with a second resonance energy fixed at the twice as a foundamental energy, and with the width of 
a second resonance fixed at the twice as a foundamental one. ClxC2 represents foundamental and second CRSFs 
with both resonance enrgy free, but with the width of a second resonance fixed at the twice as a foundamental one. 
t Referring to equation (3), and defined at 10 keV in unit of counts keV~i cm~^ s~^. 
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Table 4. Best-fit parameters of the phase-resolved 2005-1-2006 spectra of Her X-1. 
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0. / 


9 9+1.4 
— 1.2 


'•^-0.3 


on n+0.5 


^■^-0.1 


^^■"-1.3 






n R1 ("^n^i 


MP^ m 9 




9 1 
Z. 1 


1 1+0.9 
-•^■^-0.9 


in 1+4-3 




9 9+0.3 
^■^-0.1 


id.»_2.0 




1 4+1-5 


n 4Q ( A(!i\ 


(p =i.U-l.i 






















MP V PI 








7 


■^7 4 


1 Q 


S 4 

O .4: 






1 fi7 I'^n'i 


NP V PI 9 


n 7 


9 

z.o 




1 n "i+i-s 

iU.d_l g 


■^"■0-0.5 


9 q+0.2 
^■•J-0.4 


11 9+1-° 
-"^-■^•^0.9 




9 4+0-7 

^•*-l.l 


1 n7 l'4Q'l 


NPxClxC2 


0.4 


3.0 


1 7+14.4 
-■^•'-1.7 


o q+3-3 


36.5t°;^ 


9 9+0.4 
^■^-0.2 


10.21^9 


70.2t^:^ 


1 fi+°-3 

^■"-0.7 


1.08 (48) 


0=1.1-1.2 






















NPxCl 


1.9 


1.7 


-^•^-1.1 


6 e+O"* 

"■"-0.3 


34.8t°J 


1 o+"-2 

^■"-0.2 


4.5+2/ 






0.62 (50) 


NPxC12 


1.9 


1.7 


1 2+1-'' 


"■"-0.3 


34.8l°:« 


^■^-0.2 


4 t;+2.6 
+ .O2 1 




0+1-1 

"■"-0.0 


0.63 (49) 



* Abbreviations and the definition of parameters are the same as in table 3. 
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